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1.0 


INTRODUCTION 


It  is  possible  to  use  a  simple  shear  spectral 
model  to  produce  some  estimates  of  the  mean-square  velocity 
difference  over  given  values  of  The  distortion  of  a 
fluid  region  (a  wake  cross  section,  say)  of  dimension  L  is 
clearly  due  to  velocity  differences  over  4?  equal  to  L 
and  smaller.  In  this  note  we  explore  the  behavior  of 
velocity  differences  over  A2-  with  a  simple  model  for  the 
shear  spectral  density  function. 

2.0  FORMULATION 


From  (1)  it  is  apparent  that 


f(B )»  a  <u($yu*(fi)>  (3) 

is  the  horizontal  kinetic  energy  wavenumber  spectral  density 
function  (defined  over  positive^).  Using  a  trigonometric 
half-angle  formula  we  can  write  (2)  as 

Oe 

<  JuftM*)- «(«)/*> S  V  .  (4) 

o 

If  we  take  the  limit  A  we  obtain 

o© 

A*1?  s  &*€(&)  (5) 

a?-*oL  J  ^ 

where  the  left  side  of  the  equality  is  the  mean-square 

rv 

vertical  shear  5  .  Equation  (5)  demonstrates  the  relation¬ 
ship  between  the  spectral  density  of  a  quantity  and  that  of 
its  gradient.  In  this  case 

(6) 

where  <tL( 6)  is  the  shear  vertical  wavenumber  spectral 
density  function.  Relationship  (6)  is  only  applicable  in 
the  limit  A?*0.  In  general 

<t>  (0;A?)  -  ?  42'^  ( t  PA*)  (7) 

where  ^^)A2)ls  the  spectral  density  of  shear  defined  as 
the  mean-square  velocity  differences  over  42  ,  divided 

by  A iV. 


Gargett  et^  al_.  (1980)  have  assembled  a  shear 
spectrum  which  we  model  as  shown  in  Figure  1.  The 

spectrum  is  flat  to  about  10  cpm  where  it  takes  a  -1  slope 
to  about  1  cpm.  A  high- wavenumber  turbulence  dissipation 
range  exists  beyond  that  point.  We  do  not  include  contribu¬ 
tions  from  that  range  in  this  analysis.  We  represent  the 
spectrum  as 

and 

,  (8) 

-  o  , 

where 

<3#-  P.Z7  A/0_/  (io  cpm) 

and  . 

6C  =  xr  ^  **  (1  cpm) 

We  substitute  (6)  and  (8)  into  (4)  and  define  to 

obtain 


(9a) 


where 


1 


and  A3  and  =  £  g  .  We  integrate  Ij  and  I2  to 

obtain 

ini'  (t*.  **,-»)  4  __ 


(10a) 


4  [  ¥*  (***-  *Vc~  0  * 

~  )^*  (co^a^- /)  +2yp  a**.  Ayt'J 
-CiUi*) +Ci<*9 c)  , 


(10b) 


where 


<0 


(10c) 


and 


CYfcO 


-1 


OO 


( lOd ) 


are  the  sine  and  cosine  integrals  which  are  available  in 
tabular  form. 


3.0 


RESULTS 


Given  a  value  of  equations  (9)  and  (10)  yield 
an  estimate  of  the  mean-square  velocity  difference  as  a 
function  of  the  spectral  amplitude  .  We  have  computed  a 
number  of  estimates  of  which  are  shown  in 
Table  1,  and  are  plotted  in  Figure  2.  A  linear  trend  is 
evident  for  greater  than  about  5m.  The  mean-square 
velocity  difference  over  16m  is  about  12  times  the  differ¬ 
ences  over  2m;  the  root-mean-square  difference  over  16m 
(which  might  be  used  to  characterize  the  distortion  rate) 
is  about  3.5  times  the  values  over  2m.  Table  1  also  shows 


Table  1 


Mean-square  velocity  differences  versus  spacing 
Afc.  Entry  with  subscript  <3^  refers  to  mean- 
square  difference  due  to  contributions  of  com¬ 
ponents  with  wavelengths  larger  than  2.A2  .  I*  is 
the  ratio  of  this  mean-square  difference  to  the 
mean-square  difference  with  all  wavenumbers 
contributing . 


2 

it;'  zni1-> 

2.14 

$«)"<'!»!!  J*> 

1.072 

i&AOtaiVs  r 

1.00  p  0.93 

4 

5.47 

1.367 

1.20 

0.88 

5 

7.10 

1.419 

1 .22 

0.86 

8 

11.78 

1.473 

1.22 

0.83 

10 

14.90 

1.490 

1.22 

0.82 

16 

24.24 

1.515 

1 .22 

0.80 

-5- 


±&Ai)  Zl*i«l\>  which  clearly  approaches  an  asymptotic 

value  of  about  1.52.  Note  that  4-?  for  large 

42.  This  implies  that  o(.  A4”*  ;  we  observe  this  sort 

of  behavior  in  the  YVETTE  shear  estimates. 

In  order  to  characterize  the  contribution  of 
components  with  vertical  wavenumbers  less  than  any  given 
value  ft)  to  the  mean-square  velocity  difference,  (which  we 
will  denote  as  we  simply  repeat  the  above 

analysis  but  with  the  upper  limit  of  integration  in  (4)  set 

to  p  .  Of  special  interest  might  be  the  contribu¬ 

tion  to  the  mean-square  velocity  difference  over  due  to 
components  with  vertical  scales  greater  than  A3-.  We  define 
0=  77  /ta  which  means  that  we  are  considering  the  contribut¬ 
ions  to  the  mean-square  velocity  difference  o’  components 
with  wavelengths  s  *LA2  and  longer. 

The  results  of  these  calculations  are  also 

shown  in  the  table  (column  4).  The  value  of 
Z&o**)  approaches  an  asymptotic  value  of 

1.22.  Also  shown  in  the  table  is  the  ratio 

r  =  <  /zm s»*> 

The  results  show  that  93%  of  the  mean-square  velocity 
difference  over  2m  is  due  to  components  with  wavelengths 
of  4m  and  longer.  This  value  drops  to  asymptotically 
about  80%. 


4.0 


DISCUSSION  AND  A  SIMPLE  SIMULATION 


These  results,  which  are  based  on  a  spectral- 
density-function  model,  suggest  that  the  mean-square  velo¬ 
city  difference  over  values  of  AE  in  our  range  of  interest 
increase  roughly  as  AE  and  that  Sl(A3)  H.  <1(AU/1) 
behaves  as  AZ'(  .  Trends  of  this  sort  are  evident  in  YVETTE 
data.  If  we  characterize  the  relative  distortion  rate 
between  locations  A2  apart  as  ,  we  see  that  the 

distortion  rate  behaves  as  42^-  This  implies  that  the 
treatment  of  velocity  variation  as  a  linear  variation  over 
scales  of  interest,  for  wake  model  application  say,  may  be 
inappropriate.  However,  90%  or  better  of  the  root-mean- 
square  distortion  rate  of  two  points  separated  by  42  is 
given  by  components  with  wavelength;  AA2  and  longer.  This 
means  that  while  we  can  adequately  characterize  the  distor¬ 
tion  rate  between  two  points  separated  by  42  with  velocity 
data  resolved  to  wavelengths  of  2AE  and  longer,  we  should 
not  characterize  the  distortion  between  these  locations  with 
a  linear  variation. 

In  an  attempt  to  make  these  results  more  vivid,  we 
undertook  an  elementary  kinematic  "wake"  simulation  using 
the  velocity  field  recorded  during  YVETTE  deployment  09  in 
the  Western  Sargasso  (Lambert  et  aT . ,  1980).  The  velocity 
profile  is  shown  in  Figure  3;  the  data  have  an  effective 
resolution  of  about  3m.  We  computed  the  distortion  of 
initially  nearly  circular  (10m  diameter)  patterns  of  passive 
tracers  located  in  depth  intervals  70-80m,  110-120m,  220- 

230m,  345-355m,  with  the  velocity  field  fixed  over  the 

simulation  time  interval.  The  resulting  "wake"  cross- 
sections  to  one-hour  time-late  are  shown  in  Figure  4. 


> 


In  three  of  the  four  cases,  there  clearly  is  substantial 
distortion  over  scales  less  than  the  initial  wake  cross 
section . 


Finally,  we  note  that  the  work  of  Patterson  et  al . 
(1981)  with  YVETTE  data  suggests  that  **'►0  sc  , 
where  /V*  is  the  mean  (depth  averaged)  Brunt-Vaisala  fre¬ 
quency.  If  the  spectral  shear  representation  is  appropriate 
we  can  estimate  S  over  other  £?  via  (9)  and  the  estimate  of 
<po  given  * 


> 


I 
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Figure  3  YVETTE  09  Velocity  Profile  used  for  Simulations 


DISTANCE  ai> 


Figure  4  Distortion  of  nearly  circular  pattern  of  passive  tracers  in  YVETTE 
09  velocity  field.  Note  the  discontinuous  vertical  scale. 

(a)  Initial  Pattern 


at  tern  after  10  minutes 


YUETTE  09 


DISTANCE 


DISTANCE 


YUETTE  09 
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